Chaotic and Arnold stripes in weakly chaotic Hamiltonian systems 
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The dynamics in weakly chaotic Hamiltonian systems strongly depends on initial conditions and 
little can be affirmed about generic behaviors. Using two distinct Hamiltonian systems, namely one 
particle in an open rectangular billiard and four particles globally coupled on a discrete lattice, we 
show that in these models the transition from integrable motion to weak chaos emerges via chaotic 
stripes as the nonlinear parameter is increased. The stripes represent intervals of initial conditions 
which generate chaotic trajectories and increase with the nonlinear parameter of the system. In 
the billiard case the initial conditions are the injection angles. For higher-dimensional systems and 
small nonlinearities the chaotic stripes are the initial condition inside which Arnold diffusion occurs. 
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In weakly chaotic Hamiltonian systems the domains of 
regular dynamics in phase space are very large compared 
to the domains of chaotic dynamics [1] . The general char- 
acterization of the full dynamics in such systems is ex- 
tremely difficult. Starting from integrable Hamiltonian 
systems, we show that for weak perturbations the chaotic 
trajectories are solely generate inside specific intervals 
of initial conditions. These intervals increase with the 
magnitude of the perturbation creating stripes of irreg- 
ular motion, which are designed here as chaotic stripes. 
Chaotic motion for the whole system is only obtained for 
larger perturbations via the overlap of the chaotic stripes. 
Results are shown for one particle inside an open rect- 
angular billiard with rounded corner and four particles 
coupled globally on a discrete lattice. Thus, we show 
that the way weak chaos is born in the very complex 
dynamics of the coupled maps, is identical to what oc- 
curs in the simple billiard system. For higher-dimensional 
systems the chaotic stripes contain the chaotic channels 
where Arnold diffusion occurs and are designed as Arnold 
stripes. 



I. INTRODUCTION 

The precise description of weak chaos in noninte- 
grable high-dimensional Hamiltonian system is not triv- 
ial. When starting from an integrable system, the large 
amount of regular structures (invariant tori in two di- 
mensions) which remain for small nonlinear perturba- 
tions strongly influences the dynamics [2]. The phase 
space dynamics is mainly regular with few chaotic trajec- 
tories. In distinction to what happens in totally chaotic 
or ergodic systems, the dynamics is now strongly depen- 
dent, in a very complex way, on the initial conditions 
(ICs). The present contribution shows a simple and clar- 
ifying way to analyze the whole dynamics when weak 
chaos is dominant. The dynamics dependence on initial 
conditions is shown in the plot: ICs versus the nonlin- 
ear parameter. Chaotic trajectories arc generated inside 



intervals of ICs which get larger as the nonlinear parame- 
ter increases, forming chaotic "stripes" -like structures in 
the plot ICs versus nonlinear parameter. Thus the pur- 
pose of the present work is to show the way {via chaotic 
stripes) integrable Hamiltonian systems are transformed 
in weakly chaotic systems. Since chaotic stripes are ex- 
pected to exist for a large class of Hamiltonian systems, 
results are shown here for two entirely different systems. 

The first system is the two dimensional open rectan- 
gular billiard where the corners have been modified, they 
are rounded. Such open billiard, but with straight cor- 
ners and a rounded open channel, was studied recently 
[3]. It was shown that the rounded open channel gener- 
ates a chaotic dynamics inside the stripes and thus inside 
the billiard. The chaotic motion is characterized based 
on exponential decays of escape times (ETs) statistics [4] 
inside the stripes. We would like to mention some exam- 
ples of two-dimensional billiards where boundaries have 
been modified: edge roughness in quantum dots [5], soft 
billiards [6-9], edge corrections [10] in a resonator, effects 
of soft walls [11], rounded edges [12, 13], deformation of 
dielectric cavities [14], location of the hole [15], among 
others. 

The second example considered here comes from the 
nonlinear dynamics, namely 4 particles which are globally 
coupled on a discrete lattice. This system was extensively 
studied in [16-18], and the dynamics characterized by the 
existence of an ordering process called clustering. In the 
context of the present work, this model serves to show 
that chaotic stripes exist in a higher-dimensional system 
and arc related to the Arnold web, or stochastic web [1]. 
Different from the billiard case, the chaotic motion will 
be characterized by the Finite Time Lyapunov Exponent 
(FTLE). Starting from a chaotic trajectory, regular struc- 
tures induce sticky motion [19-21], the convergence of the 
FTLEs is affected and the dynamics becomes strongly 
dependent on ICs. Recent methods [11, 22-26], which 
use higher-order cummulants of Gaussian-like distribu- 
tions of the FTLEs to characterize the whole dynamics 
via sticky motion, do not work quite well for very small 
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nonlinear perturbations since the distributions may be 
multimodal [27]. 

Although we are considering two totally distinct phys- 
ical systems, the common feature is that the appearance 
of chaotic trajectories for weak perturbations occurs via 
the same process, the chaotic stripes. In other words 
we are interested to understand how the transition from 
integrable motion to weakly chaotic regime happens as 
the parameter of nonlinearity is increased. This occurs 
for the very simple billiard system and for the complex 
coupled maps, for which chaotic stripes are shown to be 
related to the Arnold web. Since the second model con- 
tains coupled standard maps, results should be valid for 
a large class of dynamical systems. 

The paper is organized as follows. Section II presents 
the model of the open billiard system and shows numeri- 
cal results for the ETs decays and the existence of stripes 
in emission angles. Section III presents the coupled map 
lattice model, shows the existence of the stripes and their 
relation to the Arnold web. In Section IV we present our 
final remarks. 



II. STRIPES IN OPEN BILLIARDS 

The model used is shown in Fig. 1, which is a rectan- 
gular billiard with dimension L x D. The four corners 
are assumed to be rounded with radius R. For R/L = 





FIG. 1. (Color online) The rectangular billiard, with dimen- 
sion Lx D, has four rounded corners (with radius R, see mag- 
nification of just one border). The escape point with aperture 
a lies exactly in the middle of the billiard. Initial angle 9o and, 
schematically, the shortest escape trajectories {P, E, A2 and 
A4) are shown. In all simulations we use L = 4 and D = 10. 

some orbits {P, E, A2, A4) are shown schematically in 
Fig. 1, whose inital angles are 6'g"'' = arctan [2^]- For 
the closed billiard these are periodic orbits with period 
g = 2" (n = 1,2,3,...). 

The closed rectangular billiard is integrable in the ab- 
sence of corner effects (i? = 0). All Lyapunov exponents 
are zero and the dynamics can be described in terms of 
invariant tori. Such regular dynamics is shown in the 
the phase-space of Fig. 2(a), where the collision angle 9 
is plotted as a function of horizontal coordinate x. Tori 
with irrational winding numbers are the straight lines 
parallel to the x-axis. Tori with rational winding num- 



bers are periodic orbits and are the marginally unsta- 
ble periodic orbits (MUFOs) from this problem. When 
opening up the billiard, ICs which start on an irrational 
torus will certainly escape the billiard after some time, 
but trajectories which start exactly on a rational tori 
will never leave the billiard, except for those ICs which 
match the opening channel. Figure 2(b) shows the phase- 
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FIG. 2. (Color online) Phase-space dynamics for (a) R = 
and (b) R/L = 0.1 (one IC): Green (light gray) points 
occur when the trajectory collides with the rounded corner 
and defines the scattering region (SR) or chaotic region. Blue 
(dark gray) points occur when the trajectory collides with the 
vertical and horizontal parallel walls. 

space dynamics when R/L = 0.1 starting from one IC 
6*0 = 0.86. Blue points are related to points along the 
trajectory which collide with vertical and horizontal par- 
allel straight walls, while green points are those which 
collide with the rounded corner. What happens is that 
the trajectory starts on an irrational torus, travels along 
this torus (blue horizontal lines) until it collides with one 
of the rounded corners (green points) and is scattered 
there to a new torus. It travels on this new tori until 
it is scattered again to another tori. This procedure re- 
peats itself until the whole phase space is filled. This is 
exactly what happens in Fig. 2(b), obtained using just 
one initial condition. The rounded corner is called scat- 
tering region (SR) and generates the chaotic dynamics. 
Although there are no regular islands in phase space, 
sticky motion occurs due to the existence of MUPOS, as 
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explained in details in [3]. For some of the horizontal 
lines in Fig. 2(b), there are some regions without points. 
These are trajectories very close to the periodic orbits 
(closed billiard case) shown in Fig. 1, and before they 
can visit all the points on the horizontal line, they are 
scattered by the SR to another horizontal line. 

In order to analyze the effect of rounded corners on the 
emission angles dynamics we start simulations at times 
t ~ from the open channel, uniformly distributed from 
xo = D/2, with an initial angle 0o towards the inner part 
of the billiard with velocity \v\ = 1. Elastic collisions are 
assumed at the billiard boundaries and at the rounded 
corners of the open channel. For each of the 10^ IC dis- 
tributed uniformly in the interval 1.00 x 10~^ ^ ^ 1-57 
we wait until the particle leaves the billiard and record 
6f. In all cases a = 0.4. 

The complex dynamics generated by the rounded cor- 
ners becomes apparent when the escape angles Of are 
plotted as function of the initial incoming angle 6q and 
for different ratios R/L. This is shown in Fig. 3 which 
was obtained by using 10^ x 10^ points in the interval 
2.0 X 10-3 < R/L < 0.35 [-6.0 < log(i?/L) < -1.0]. 
Figure 3(a) shows 9f as a function of log(i?/L) and 9q. 
Each color is related to one emission angle 6 f (see the col- 
orbar on the right: dark blue — >■ red — yellow — > white in 
the colorscale and dark gray to white in the grayscale). 
These emission angles vary between 6*/ ^ 1.4 (almost 
horizontally to the left) and 0/ ~ 4.5 (almost horizon- 
tally to the right). It is possible to observe [better seen 
in the magnifications. Figs. 3(b) and (c)] that there are 
always intervals of ICs which leave to the same Of (same 
color). As R/L increases, they appear as horizontal lines 
with the same color, and are called here as "isoemis- 
sions stripes". As R/L increases more and more, some 
isoemissions stripes survive while others are destroyed or 
mixed. Examples of such stripes are seen in Fig. 3(a) 
around ^ 0.7 (yellow), in Fig. 3(b) around Oq ^ 0.9 
(brown), among many others. In general, the emission 
angles show a very rich dynamics due to the increasing 
rounded corners, alternating between all possible colors. 

Figure 3(b) shows a magnification from Fig. 3(a). The 
magnification is taken around the isoemissions stripe re- 
lated to the trajectory E {Oq ^ 0.89) shown in Fig. 1. We 
observe that a new stripe emerges symmetrically from 
the middle of the isoemissions stripe. Different from the 
isoemissions stripe, inside this new stripe inumerous dis- 
tinct escape angles (colors) appear. This new stripe is 
called "chaotic stripe" for reasons which become clear 
later. The width of the chaotic stripe increases linearly 
with R/L. Outside the isoemissions stripe a large amount 
of smaller stripes with distinct angles can be observed. 
In order to show this in more details and to explain the 
physics involved in the emission angles, we discuss next a 
magnification of Fig. 3(b) (see white box). This is shown 
in Fig. 3(c), where a sequence of isoemissions (yellow and 
purple) and chaotic (all colors) stripes appear. ICs which 
start inside the chaotic stripes collide, at least once, with 
the rounded corner so that distinct emission angles can 
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FIG. 3. (Color online) (a) Escape angle as a function of 
log (-R/L) and 9o, (b) and (c) are magnifications. 



be observed. Different from what is observed from bor- 
der effects [3], here the sequence of isoemissions stripes 
in Fig. 3(b), and the corresponding multicolor chaotic 
stripes from Fig. 3(c), are not born at the boundary 
between the purple (black) and orange (white) escape 
angles at R/L ~ 0. The rich and complex dynamics 
is always generated from ICs starting inside the chaotic 
stripes, so that a large amount of different colors in the 
emission angles appears, showing that tiny changes or er- 
rors in the initial angle may drastically change the emis- 
sion angle. The location of the chaotic stripes itself is not 
self-similar, but inside the chaotic stripes the self-similar 
structure is evident. Chaotic stripes which emerge at 
different initial angles at R/L ^ 0.0 increase their width 
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with R/L and start to overlap for R/L ~ 0.05, where the 
dynamics of strongly sensitive to initial angles ^o- The 
width of the stripes increase with R/L and allows us to 
characterize the different domains of the dynamics. 

The dynamics generated by starting with IC inside the 
chaotic stripes is the consequence of trajectories which 
collide with the rounded corner and the chaotic motion 
appears. The collision with the rounded border was 
checked numerically (not shown) and the chaotic motion 
has to be demonstrated. To do this we use the ETs statis- 
tic defined [28] by Q{t) = limjv-i-oo where A'' is the 
total number of recurrences and N^- is the number of re- 
currences with time T > t. The time T is recorded when 
the trajectory returns to the recurrence region, which is 
the open channel. In our model, when a trajectory re- 
turns to the open channel it leaves the billiard. Thus, 
since we have an open system, the number of recurrences 
is counted over distinct ICs. In chaotic hyperbolic sys- 
tems the ETs statistic decays exponentially while systems 
with stickiness it decays as a power law Q{t) oc r~'''=='= 
with 7esc > 1-0 being the scaling exponent. Usually it 
is assumed to have stickiness when a power-law decay 
with 7esc > 1.0 is observed for two decades in time. 
Power-law decays with jesc = 1-0 [3, 29] appear in in- 
tegrable systems, which is the case of our model when 
R/L = 0.0. Figure 4 displays in a semi-log plot the 
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FIG. 4. (Color online) Semi-log plot of Q(t) for distinct values 
of R/L. Each curve was obtained using 2.0 x 10'' ICs. 

quantity Q{t) obtained using different values of the ratio 
R/L = 0.0, 1.0 X lO"'^, 1.0 X 10-3, 1.0 X lO^^, 5.0 x 10-^, 
where L is kept fixed. ICs arc 1.00x10"^ < Oo < 1.57 and 
-1.00 x 10"^ < 6*0 < -1.57. Straight lines in this plot are 
exponential decays. For the integrablc case R/L = 0.0 
(curve e) a power law decays occurs with jegc ^ 1-0. Tra- 
jectories close to the periodic orbits inside the rectangle 
generate the power-law decay with j^sc = 1-0. Small cor- 
ner effects, R/L = 1.0 x 10~* 1.0 x lO""^ (see curves 
d — >■ c), change the qualitative behavior of Q{t) so that 
regular trajectories close to MUFOs generate, for small 
period of times, a power-law decay with 1.0 < 'jesc < 2.0. 
For a detailed discussion of a similar behavior we refer the 



reader to [3]. However, for R/L = 1.0 x 10-2[log(i?/L) = 
-4.6], 5.0 X 10-2[log(i?/L) = -3.0] (curves b,a), which 
correspond to values from Fig. 3(c), a straight line decay 
is observed, characterizing the transition to an almost 
totally chaotic motion when stripes overlap. 

Summarizing results for the billiard model, we showed 
that specific intervals of ICs, called chaotic stripes, which 
increase linearly with the nonlinear parameter, generate 
chaotic trajectories. Outside these stripes we have the 
isoemissions stripes, for which the ICs generate a regular 
motion. For smaller values of the parameter, we have 
a mixture of regular and chaotic motion, i.e. a mixture 
of a large amount of isoemissions stripes and few chaotic 
stripes, and the system is weakly chaotic. Only when the 
nonlinear parameter increases, so that all chaotic stripes 
overlap, the whole system becomes chaotic. 



III. STRIPES IN COUPLED MAPS 

At next we analyze the IC dependence of the dynamics 
of 4 particles globally coupled on a discrete lattice. It has 
a 8-dimensional phase space and particles are coupled on 
a unit circle, where the state of each particle is defined 
by its position 27ra;'-*' and its conjugate momentum p^'^\ 
The dynamics on the lattice obeys [16] 
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where mod 1 must be taken in the variables p^'^ , a;^*^ , and 
i = 1, . . . ,N = 4. For K > the interaction between 
two particles i and j is attractive [17]. For this coupling 
the total momentum Px = '^JSi^Pt is preserved. This 
model has 8 Lyapunov exponents which come in pairs. 
The reasons for using such a model are: (i) these are 
coupled standard maps, thus their properties should be 
valid for a large class of physical systems; (ii) the phase 
space dynamics is 8, thus it can be used to show that, 
in case results are similar to those from Section II, they 
are valid for higher dimensions and connect the chaotic 
stripes to the Arnold web. 

An interesting way to describe the complicated dynam- 

(2) 

ics dependence on ICs is shown in the plot K x Pq from 
Fig. 5. In colors is plotted the largest FTLE after 10^ 
iterations. Initial conditions are chosen on the invariant 
structure Pt = 0.0 and Xcm ~ 0.0. This plot shows, as 
a function of K, those ICs Pq^'' which generate the regu- 
lar, mixed or chaotic motion. While black colors describe 
initial conditions with generate a regular motion, with 
zero FTLEs, cyan, red to yellow colors are obtained from 
ICs which generate increasing positive FTLEs. A rich 
variety of structures is observed. Essentially two distinct 
larger regions of motions are observed and demarked in 
the plot as: (A) where FTLEs are smaller and (B) where 
FTLEs arc larger. In between the FTLEs mix themselves 
along complicated and apparently fractal structures. For 
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FIG. 5. (Color online) Grid of 1000 x 1000 points with the 
FTLEs as a function of K x p^^^ . 



random states, the regular structures affect locally the 
FTLEs inducmg sticky motion, so that the correspond- 
ing FTLEs decrease. Thus each point in the mixed plot 
from Fig. 5(a) with and intermediate FTLE for a fixed K 
value, is necessarily related to a trajectory which suffers 
sticky motion. Figure 5(b) is a magnification of Fig. 5(a) 
and shows similar stripes structures observed in the bil- 
liard case. As the nonlinear parameter K increases, the 
width of the chaotic stripes increase and they start to 
overlap, generating the totally chaotic motion, but only 
for higher K values (not shown). 

It is very interesting to observed that the chaotic 
stripes define the channels of chaotic trajectories which 
form the Arnold web [1] which allows the Arnold diffu- 
sion to occur in higher-dimensional systems. This can 
be better understood when the FTLEs arc plotted in the 
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= 0.15 [see dashed line in Fig. 5 (a)] 



phase space projection ' x Pg^'', shown in Fig. 6 for 
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if — >■ FTLEs also go to zero inside region (A). How- 
ever there are some stripes, which emanate from K ^ 
for which FTLEs arc larger. Two larger stripes are born 

(2) 

close to Pq 0.0 and growth symmetrically around this 
point as K increases. ICs which start inside the stripes 
generate chaotic trajectories. Thus these stripes corre- 
spond to the chaotic stripes found in the billiard case in 

(2) 

Sec. II. Along the line Pq ^ 0.0 the dynamics is almost 
regular for K < 1, becoming chaotic for larger values 
of K. This picture shows us clearly that for the same 

(2) 

K value, different ICs (in this case pg ) have distinct 
FTLEs. Zero FTLEs are related to regular trajectories, 
larger FTLEs are related to chaotic trajectories and inter- 
mediate FTLEs are related to chaotic trajectories which 
touched, for a finite time, the regular structures from the 
high-dimensional phase space. Such regular structures 
can be global invariants (as Pt) or local, or collective 
ordered states which live in the high-dimensional phase 
space. As trajectories itinerate between ordered and 



FIG. 6. (Color online) FTLEs in the phase space projection 
p^^' xpj,^' for K = 0.15 [see white dashed line in Fig. 5(a)]. 

are related to close to zero FTLEs, while purple, red to 
yellow points are related to increasing FTLEs. This pic- 
ture is similar to the one shown in [18] (but in black and 
white), and presents a form called as "onion". In fact, 
this plot shows an Arnold web, or stochastic web, where 
points with larger FTLEs are the chaotic channels inside 
which the Arnold diffusion occurs. The FTLEs along the 
line p'lP — 0.0 in Fig. 6 are exactly the FTLEs along 
the white dashed line in Fig. 5(a). In other words, ICs 
for K ~ 0.15 [white dashed line in Fig. 5(a)], inside the 
chaotic stripes, are those which generate trajectories with 
positive FTLEs and correspond exactly to the positive 
FTLEs from Fig. 6 for pg^' = 0.0 (see also white dashed 
line). Thus the chaotic stripes define the width of the 
chaotic channels, for a given isT, inside which Arnold dif- 
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fusion occurs as time increases. Summarizing, the chaotic 
stripes in Fig. 5 show directly the intervals of ICs where 
the chaotic channels are located and how the Arnold web 
changes with K. 

IV. CONCLUSIONS 

Chaotic stripes are shown to generate weakly chaotic 
dynamics for two distinct Hamiltonians systems, namely 
an open rectangular billiard with rounded corners and 
4 globally coupled particles on a discrete lattice. The 
width of the stripes increase with the nonlinear param- 
eter and allow us to characterize the different domains 
of the dynamics. For the billiard case the stripes define 
injection angles which generate an chaotic dynamics in- 
side the billiard and an uncertainty about ejection angles. 
For small values of the nonlinear parameter but higher- 
dimensional systems, the stripes are the chaotic channels 



in the Arnold web inside which the Arnold diffusion oc- 
curs [1] . From this perspective we can define the chaotic 
stripes as Arnold stripes. We also observed such stripes 
for Hamiltonian systems with higher values of > 4 and 
other couplings (not shown here). This suggests that 
weak chaos in Hamiltonian systems always emerges in- 
side stripes of particular ICs. This is of relevance since 
in such systems finite time distributions, independently 
of the considered physical quantity, are not Gaussians 
anymore and decay rates of recurrences and time corre- 
lations are difficult to obtain [24, 25]. For higher values 
of the nonlinear parameters a totally chaotic motion is 
observed due to the overlap of Arnold stripes. 
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